Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Metallic State of the Three-band Hubbard Model with Super-lattice 

Structure 

Shigeru Koikegami^ * , Takashi Yanagisawa^ , and Masaru Kato^ 

^ Nanoelectronics Research Institute, AIST Tsukuha Central 2, Tsukuha 305-8568, Japan 
'^Department of Mathematical Sciences, Osaka Prefecture University, Sakai 599-8531, Japan 

(Received January 27, 2006) 

We investigate the dynamical superlattice correlation in the two-dimensional three-band 
Hubbard model on the basis of the unrestricted fluctuation exchange approximation. We 
calculate the one-particle spectral function, the spin correlation function and the charge 
correlation function at finite temperature. We find that some experimental results can be 
reproduced consistenly by taking inhomogenous distribution of Cu 3d electrons into account. 
The correlation functions suggest that several kinds of instabilities with spatial inhonioge- 
nieties exist in some regions, where these instabilities significantly affect the one-particle 
spectral fuctions. 
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correlation effect, spatial inhomogeneity 

1. Introduction 

The quasi-one-dimensional (QID) charge order, which is well-known as the stripe order, 
has been one of the significant issues in high-Tc cuprates (HTC) for recent years. The incom- 
mensurate peaks observed by the neutron scattering in La2-xBaxCu04 are to be explained 
on the assumption that the static stripe order is present around 1/8-filling.^ The presence of 
the static stripe order in the ground state has been proved theoretically by the density ma- 
trix renormalization group method on the basis of the QID t-J model. ^ Originally, since the 
beginning of HTC studies many experimentalists have pointed out the so-called 1 / 8-problem 
in La2-xBaxCu04, i.e., the remarkable suppression of superconducting transition tempera- 
ture at X ~ 1/8. Many theorists have thought that the strong on-site Coulomb interaction 
makes electrons almost localized with the long-period spin and charge correlations around 
1/8-filling. They have studied these long-period correlations by the numerical or analytical 
methods on the basis of the two-dimensional (2D) Hubbard models and have clarified that 
some long-period correlated states can be the ground states around certain fillings. ^"^■^ These 
theoretically predicted long-period correlation states are consistent with the stripe state de- 
duced from the neutron scattering experiments.^^ 

The theoretical results on the basis of the 2D Hubbard models mean that the strong 
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on-site interaction is enough to stabilize the stripe state near l/S-fUhng as far as the ground 
state is considered. However, at finite temperature the strong fluctuations characteristic to 
the low dimensional strongly correlated electron systems exist. These strong fluctuations are 
composed of the many elements, e.g., the anti-ferromagnetic (AF) fluctuation and the charge 
ordering fluctuations with certain classes of super-lattice structures. These strong fluctuations 
accompany with the instabilities of the spin and charge orderings and simultaneously tend to 
interrupt them. The stripe with such strong fluctuations in metallic state is called a dynamical 
stripe. ^^"^'^ The angle resolved photo-emission spectroscopy (ARPES) shows that the stripes 
in Bi2Sr2CaCu208+5 has more dynamic character than those in La2-xSrxCu04. The discrim- 
ination whether the stripe is static or dynamic is to be judged from the characteristic energy 
of the fluctuation mode. If the fluctuation mode has low characteristic energy, i.e., the stripe 
fluctuates slowly, the mode can easily couple with phonon and the static stripe long-range 
order will appear. On the other hand, if the fluctuation mode has high characteristic energy, 

1. e., the stripe fluctuates fast, the mode cannot couple with phonon and the system will stay 
metallic together with the fluctuation. When we consider these fluctuations in the analysis of 
these models without the long-ranged Coulomb interaction, we can expect that the metallic 
state should be recovered at higher temperature. In this state, however, the electronic prop- 
erty is to be still affected by the stripe instability as far as its dynamical factor is investigated 

In this paper, we analyze the 2D three-band Hubbard model on the basis of the unre- 
stricted fluctuation exchange approximation (UFLEX). UFLEX can consider a strong corre- 
lation effect and a number of the spatial inhomogencities without any assTimptions aboT.it the 
ground state. These considered effects contain not only a static component but a dynamic 
one. As a result, the electronic state reproduced on the basis of UFLEX is affected by a num- 
ber of spatially inhomogeneous fluctuations simultaneously. In fact, when we calculate the 
one-particle spectral functions and the spin correlation functions for our fully self-consistent 
solutions, we can find that the one-particle spectra of the electronic state reflect these exotic 
state accompanied with the anti-ferromagnetic or the further long-periodic instabilities. Es- 
pecially, near 1/8-filling the spin and charge correlation functions show the existence of the 
QID fluctuation at finite frequency. This fast QID fluctuation is originated from the striped 
spatial inhomogeneity. Thus, our results suggest that the spatial inhomogeneities should be 
concerned in the strongly correlated electron system when the electronic behavior at a finite 
frequency is discussed. 

2. 2D three-band Hubbard model 

Our 2D three-band Hubbard model Hamiltonian, H, is composed of d-electrons at each 
Cu site and p-electrons at O site. We consider only the on-site Coulomb repulsion U between 
(/-electrons at each Cu site. Then, H is divided into the non-interacting part, Hq^ and the 
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interacting part, Hi, as 

H = Ho + Hi - /.^(4.rfk. +fklfL)- (1) 

Here dko-('^ko-) ^^"^ PtLo'\p\}a^^ ) annihilation (creation) operator for p- and 

electron of momentum k and spin a, respectively, is the chemical potential. The non- 
interacting part Hq is represented by 



pL 



(2) 



V -CI Cl J 

where A^p is the hybridization gap energy between d- and p-orbitals. We take the lattice 
constant of the square lattice formed of Cu sites as the unit of length, and we can represent 
^x{y) _ 2H^^ '^^M. and Ck = — 4tpp sin ^ sin ^ , where t^p is the transfer energy between a 
d-orbital and a neighboring p^^^/) .orbital and tpp is that between a ;?^-orbital and a p^-orbital. 
In this study, we take t^p as the unit of energy. The residual part, Hi, is described as 

^1 = 4+qT4'-qi'^k'iC?kT- (3) 

kk' q 

where U is the on-site Coulomb repulsion between d-orbitals and N is the number of k-space 
lattice points in the first Brillouin zone (FEZ). 

3. UFLEX 

We have developed UFLEX in order to analyze the spatially inhomogcncous system. In our 
UFLEX formulation wc introduce 'cluster' momentum besides 'usual' momentum. The notion 
of this cluster momentum is suggested from the one in dynamical cluster approximation, in 
which the nonlocal correlations are considered beyond dynamical mean field theory. ^^"^"^ 

However, our cluster momentum is not directly related with a coarse graining of the 
Brillouin zone. When the most stable state of our system is spatially homogeneous, only a 
cluster momentum, K = 0, is enough to obtain the true solution for our problem. In that 
case, UFLEX is consistent to 'usual' FLEX.-*^^ The advantage of our introduction of cluster 
momenta is its flexibility in the correspondence to the problem. If the problem have some 
important classes of the spatial inhomogeneities, we will only need to provide the cluster 
momenta corresponding to those classes. For example, when the most important instability 
is the anti- ferromagnetic one, two cluster momenta, K = and K = {^,^), are enough. 
Of course, the stripe instabilities have their characteristic cluster momenta. Thus, when we 
anticipate to encounter the stripe instabilities in our problem, we need to provide only those 
cluster momenta. 

First, we define the unrestricted perturbed Green function, G^j^(k, ie„). Here we use the 
abbreviation of Fermion Matsubara frequencies, = 7rT(2n + 1) with integer n, where T 
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is the temperature. K indicates the cluster momentum. The set, {K}, is chosen so that any 
sum, k + K, exists in the set, {k}. GJj^(k, ie„) is calculated by the Dyson equation : 

[GjK(k,ie„)]-' = -^K ■ {G;(°\k,i6„)}"' - Sic(k,ie„), (4) 

where (5k means Kronecker's delta, i.e. (5k=o = 1 and Sj^^q = 1. [• ■ ■]~^ indicates the inverse 
operation defined so as to satisfy the identity : 

Sk = Y1 GlK-K'(k + K' - K, i6„) [GjK'(k, ien)]-' (5) 

K' 

for all k and n. G^'''^\k,ien) in Eq. (4) is the unperturbed Green function derived by diago- 
nalizing of H — Hi as 

{G;(°)(k,z-M)}"' = 

2z{2 — cos kx — cos ky) + 8tpp(l — cos kx){l — cos ky) 
^~ * z2-4t2p(l-cosA;^)(l-cos ky) ' ^ ' 

In order to estimate our unrestricted self-energy, SK(k, ie^), in Eq. (4), we adopt the UFLEX 
as follows. 

Sfc(k,ien) = 

'T'TT 

^EG<rK(k-q,i6„0e'^"'1,_^o 

qn' 

+ § E GjK-K'(q-k + K'-K,ia;^-ie„)l>P)(q,ia.^) 

q K' m 

T 



qK'm 

+ <^dK-K'(k - q, ien - iw^jF^/^^^^ (q, iw^) , (7) 



C^^'?!'K(q,iwm) 

- t^' E -^K-K' (q, ic^m) [ 5k' + (q + K - K', ia;^)] , (8) 



K' 



+ U^Y1 XKlK'(q. i^m) [<5k' - C/X^c7"(q + K - K', iu;^)] , (9) 

K' 
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and 

K' 

X [ <5k' - U^J2 Xk"-k" (q + K - K', ium) 

K" 

X XK'""(q + K-K",ia;^)]-\ (10) 
where Um = 2m 'kT with integer m are Boson Matsubara frequencies, 

^5'K(q,i'^m) = 

Grf^K-K'(q - k, io^m - ien) G^K'(k, ien), (11) 

kK'n 

and 

XK'(q,iwm) = 

E G5K-K'(q + k,ia;^ + ie„)G5K'(k-K',ie„). (12) 

kK'n 

In order to obtain the solution satisfying the conserving law, we need to solve 
Eqs. (4,7,8,9,10,11, 12) fully self-consistently. This self-consistent procedure is shown in Fig. 1 
diagrammatically. 

In numerical calculations we divide the FBZ into 16 x 16 meshes. We take 8x8 cluster 
momenta so that we can reproduce a homogeneous state, inhomogeneous states with 2-, 4-, 
or 8-lattice period charge density wave (CDW) or spin density wave (SDW) along a(x)-axis 
or 6(y)-axis, or an AF state. We prepare 2^^ = 2048 Matsubara frequencies for temperature 
T = 0.020 ~ 180K. The other parameters : tap = 1.0 ~ 0.80 eV, tpp = 0.60 ~ 0.48 eV, and 
^dp = 0-70 ~ 0.56 eV, U = 4.0 ~ 3.2 eV, which are all common for our results. 

4. Fermi surface evolution with hole- or electron-doping 

We have obtained the fully self-consistent solutions for 0.000 < 5h < 0.057, 0.112 < 
Sh < 0.231, and 0.000 < 5e < 0.249, where 5h = 1 — rid — rip and Sc = ri^ + rip — 1. 
and rip represent the number of d- and p-electrons, respectively. We could not obtain any 
convergent solution for 0.057 < 6^ < 0.112. According to the elaborated neutron scattering 
experiments on La2-x-yBaxSryCu04,^^"^^ in this region the inhomogeneous states with longer 
than 8-lattice period SDW along diagonal, which is called diagonal stripe state, could be 
realized. Our present calculation would be short of provided cluster momenta to reproduce 
the inhomogeneous states with the diagonal fluctuations. Hereafter, we restrict our discussion 
to the states with the vertical fluctuations. 
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Fig. 1. The self-consistent procedure in UFLEX. i represents the iteration count, and e is the small 
value. 



In this section we investigate the Fermi surface evolution with doping in detail. In order 
to determine the Fermi surface, we define the one-particle spectral weight as 

A-(k,S)^-llm Gco(k,i6n)I,,„^^, 



where 



= {z[z - Adp - Sg(k, z)] - 2(1 - cos k,,y)} 

x{[z^- 4tlp{l - cos/c^)(l - cosky)][z - Adp - SS(k,z)] 
—2z(2 — coskx — cosky) — 8tpp{l — cosfcj;)(l — cosfcj^)}"^ 



(13) 



(14) 



In Eq. (13) we use Padc approximation for the method of analytic continuation. In order to 
determine the Fermi surface on the basis of our calculated results, we should calculate the one 
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particle spectral weight at E = 0, A'^(k, 0), and find the k-points at which A'^(k, 0) becomes 
large. Because the thermal fluctuation with CDW, SDW, or anti-ferromagnetic instabilities 
can generate the branched energy dispersion, we show our calculated results as follows. We 
choose the both k-points on which the one-particle spectral weight at = has the largest 
and the second largest peaks for fixed kx as far as < fcj^ < kx < vr. And we repeat this 
operation for the other parts of FBZ. In the following figures, the points on which the one- 
particle spectral weight at E = has the largest and the second largest peaks are to be 
indicated with black and gray circles, respectively. Our results have been obtained at high 
temperature, and the thermal fluctuations make the Fermi level vague. However, if we have 
empty spaces surrounded by these points, we can expect the branched energy dispersion 
near the Fermi level. The branched energy dispersion is made by a spatial inhomogeneous 
instability. Furthermore, these points can suggest the Fermi surface connectivities at lower 
temperature. 

Firstly, we should mention our results for the hole-doped states around 1 /8-filling, in which 
the dynamical stripe state might be realized. As shown in Fig. 2, near 1/8-filling we have 
eight empty spots around the intersections of the planes indicated by dashed and dotted lines, 
which relate to the anti-ferromagnetic and the 4-lattice period CDW instabilities, respectively. 
It means that the one-particle spectral weights are lost due to the fluctuations both of the 
anti-ferromagnetic instability and of the 4-latticc period CDW instability. These instabilities 
occur when some parts of the Fermi surface are connected with one another by the nesting 
vectors. The parts of Fermi sinfacc connected with one another by the nesting vectors, (tt, vr) 
or (vr, — vr), arc called 'hot spots'. On hot spots the quasi-particlcs play important roles in 
the transport phenomena. Therefore, we expect some anomalous electronic behaviors near 
1/8-filling, when the one-particle spectral weights are lost on hot spots. This is the reason 
why so-called 1/8-problem occurs around 1/8-filling. 

Secondly, we show our results for lightly-hole-doped states, in which the anti-ferromagnetic 
instability might exist. In this case we have no more eight empty spots but four empty spots 
centered on (±7r/2, ±7r/2) as shown in Fig. 3. These empty spots are called 'Fermi arc's, 
which is believed to appear when the strong anti-ferromagnetic fluctuation exists. In fact, 
these Fermi arcs have been observed in Lai,97Sro.o3Cu04 by ARPES and interpreted as the 
evidence of the anti-ferromagnetic long-range order.^^ In our calculating results, the spectral 
weights in these Fermi arcs are not completely zero, and the anti-ferromagnetic long-range 
order does not exist. However, the existence of these Fermi arcs should be originated from 
the strong anti-ferromagnetic fluctuation, evolving toward the anti-ferromagnetic long-range 
order. 

When the doped holes become less, electrons should localize due to strong coulomb re- 
pulsion among them in real materials. However, in our three-band Hubbard model electrons 
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Fig. 2. Fermi surface evolutions near 1/8-filling. Left sides : (a) for 5h = 0.112. (b) for 5h = 0.127. 
Right sides : The gray-scale images interpolated from raw data are attached on account of the 
guide to the eye. The brighter spots have the larger spectral weights, (c) Illustration of nesting 
vectors. Dashed lines represent the planes connected with one another by the nesting vectors, 
(tt, tt) or (tt, — tt). Dotted lines represent the planes connected with one another by the nesting 
vectors, (7r/2,0) or (0, 7r/2). The illustration of nesting vectors is also applicable to the other 
figures representing Fermi surface evolutions. 



cannot localize even around near half-filling and the Fermi surface are remained, as shown 
in Fig. 4. If we could reproduce the Fermi surface evolution for the whole doping region, the 
one electron spectral weight at the Fermi level should be almost vanished near half-filling. 
This result means that our formulation on the basis of UFLEX is not efficient to describe 
the Mott-Hubbard metal-insulator transition. Unfortunately this is the side-effect of UFLEX, 
in which we overestimate the thermal fluctuations and recover the itinerancy of electrons. 
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Fig. 3. Fermi surface evolutions for lightly-hole-doped states, (a) for = 0.039. (b) for 6h = 0.057. 



It would be our future issue to describe the localized electrons but with strong fluctuations 
appropriately. 

Our results for electron-doped states will be highlighted. The ARPES have already clarified 
the doping dependence of the Fermi surface in Nd2-xCexCu04±5.^^^^ Except for completely 
vanished Fermi surface near half-filling, their observed evolution of Fermi surface with electron 
doping are consistent with our calculated result. As shown in Fig. 5, in lightly-electron-doped 
region the eight pockets emerge around the cross-sectional point of Fermi surface with FBZ 
boundaries. This reflects that the vertical charge fluctuation exists in the lightly-electron- 
doped region as well as near 1/8-hole-doped region. We can guess that the vertical charge 
fluctuation gets to play an important role when the localized electron recovers its itinerancy 
as doped electron increases. When the doped-electron density density increases more, each 
two of eight pockets on the same Fermi surface section come close to each other and unite into 
one as shown in Fig. 6. Thus, in heavily-electron-doped region we have four pockets around 
the cross-sectional point of Fermi surface with kx = iky. This means that the diagonal spin 
fluctuation, which is not completely anti-ferromagnetic, is strengthened instead of the vertical 
fluctuation. This diagonal spin fluctuation can cause pseudo-gap phenomena in electron-doped 
high-Tc as observed at low temperature.^^ It contrasts with the case of the hole-doped region, 
and the difference should be caused by both a large tpp. Such a tpp makes the big difference 
of the energy dispersions between the hole-doped region and the electron-doped region. If 
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Fig. 4. Fermi surface evolutions for near half-filling, (a) for 5h = 0.022. (b) for half- filling, (b) for 
= 0.021. 



the electron were localized, the difference of the energy dispersions should appear as the 
different charge distributions. In our results the electron could not be localized, but the charge 
fluctuations appear in the different way between the hole-doped region and the electron-doped 
region instead. Hence, we can insist that the inhomogeneity of the strongly correlated system 
should be considered. 

5. Spin and charge correlations 

In this section we show how the Fermi surface evolution with doping is related to the 
development of the spin and charge correlations. The momentum dependences of spin and 
charge correlation functions reflect the spin and charge spatial inhomogeneities, respectively. 

At first we sum up our obtained results about the static spin correlations as both the in- 
commensurabilities and the intensities of the peaks in the calculated elastic neutron scattering. 
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Fig. 5. Fermi surface evolutions for lightly-electron-doped states, (a) for 6e = 0.045. (b) for de = 0.066. 
(b) for (5e = 0.086. 

They are shown in Fig. 7, where 

Sid) ^ ^xt^M [<5k - C/xic"(q - K,0)]-' . (15) 

K 

In the hole-doped region the incommensurabilities 5inc, defined by 

(5inc = Max{[Qinax - (tT, Tt)]^, [Q^ax " (tT, Tr)]^^}, (16) 

are almost proportional to (5h. This has been already observed by the neutron scattering 
experiment on La2-x-yBa,xSryCu04.^^'^^'^° The ratio of (5inc to (5h in our calculating results 
are almost half of the one in these experiments. However, we believe that this deviation is 
originated from the difference between the real hole number 6^ introduced in Cu02 plane, 
which is correspond to the one in our results, and the instoicheiometric number S used in the 
analysis in Ref. 30. 
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Fig. 6. Fermi surface evolutions for heavily-electron-doped states, (a) for (5e = 0.106. (b) for (5e 
0.144. (b) for (5e = 0.200. 
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Fig. 8. The momentum dependences of -D(q) at hole-doped regions. 

On the other hand, in the electron-doped region we cannot recognize such a relationship 
between Smc and Se- It suggests that our spin correlation functions do not simply reflect the 
Fermi surface nesting but rather the existence of inhomogeneities as discussed in Section 4. 

In order to clear our insistence, we calculate the static charge correlation functions at 
some regions. The momentum dependences of these correlation functions are shown in Figs. 8 
and 9, where 

^(q) = § E GjK'(q-k,-ie„)G,^K.(k,ie„). (17) 

kK'n 

Focusing on Fig. 8(b) and Fig. 9(b), we can recognize that -D(q) are slightly enhanced around 
the symmetric lines, qx = and qy = 0. These momentum dependences of I?(q) suggest that 
the two collective modes with the momenta (±1,0) and (0, ±1) exist. These collective modes 
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Fig. 9. The momentum dependences of -D(q) at electron-doped regions. 



can be translated as the instabilities toward vertical charge ordcrings along x-axis and y-axis, 
respectively. In fact, when d\i = 0.127 and 6c = 0.066, the Fermi surfaces have eight pockets 
around (ibTr, ib7r/4) and (ib7r/4, zbvr) as shown in Fig. 2(b) and Fig. 4(b), respectively. As 
discussed in 4, these anomalous Fermi surfaces are caused from the charge inhomogeneities, 
which appear as shown in Fig. 8(b) and Fig. 9(b). 

6. Conclusion 

In this paper, we have obtained fully self-consistent solutions for the 2D three-band Hub- 
bard model on the basis of UFLEX, in which the inhomogeneous distribution of d-electrons 
are allowed. Both in hole-under-doped and electron-doped regions the one-particle spectral 
weights behave anomalously. Furthermore, the spin and charge correlation functions reflect 
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the existence of super-lattice structures around 1/8-filling. Our numerical solutions show that 
the metallic state with spatially inhomogcneous spin or charge fast fluctuation can exist at 
high-temperature. The electronic state in one of our microscopically derived solutions cor- 
responds to the dynamical stripe state mentioned in some pioneering works. In our 2D 
model, the electrons cannot have the long-range order at finite temperature because of the 
strong fluctuation characteristic of low-dimensionality. However, in three-dimensional real ma- 
terials, the spatially inhomogcneous fluctuation may tend to have the long-range order and 
form the striped state, observed in the neutron scattering experiments. This long-range order 
formation in three-dimensional real materials might not be reproduced by our simple model 
in this paper. The existences of the phonon or the long-range Coulomb potential might have 
an important role for these long-range order formation. The theoretical researches concerning 
these complicated factors are to be expected as the future problem. 
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